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We present a general group-theoretical scheme for the transport in monolayer transition-metal 
dichalcogenides. The theory elucidates the major momentum and spin relaxation processes for 
electrons, holes and hot excitons. We expatiate on spin flips induced by flexural phonons and show 
that the spin relaxation is ultrafast in free-standing membranes while being mitigated in supported 
membranes. This behavior is universal in 2D membranes that respect mirror symmetry and it leads 
to a counterintuitive inverse relation between mobility and spin relaxation. 

PACS numbers: 



Single-layer transition-metal dichalcogenides (SL- 
TMDs) put together exotic charge, spin and valley elec- 
tronic phenomena [l|, Recent advances in the exfo- 
liation and fabrication of these materials have sparked 
a wide interest in their d-band semiconducting behavior 
and spin- valley coupHng Room-temperature mo- 

bility of the order of 100 cm^/V-s in single-layer n-type 
M0S2 transistor was demonstrated and analyzed 0-0- In 
addition, the unique time-reversal relations between spin 
and valley degrees of freedom were studied from the cir- 
cular polarization of the exciton photoluminescence (PL) 
lol-flit . In spite of these recent advances and supporting 
studies of band-structure and phonons parameters [ij- 
[Tot , there is still no unifying description of the spin and 
charge transport in these intriguing materials. 

In this Letter, we present a group-theoretical analy- 
sis of SL-TMDs. The theory elucidates pertinent mo- 
mentum and spin relaxation mechanisms, and places the 
physics of charge, spin and exciton transport on a firm 
footing. We first delineate the transport limitations at el- 
evated temperatures via zeroth-order selection rules and 
make connection with the energy relaxation of hot exci- 
tons. It is shown that spin-conserving scattering between 
direct and indirect exciton bands leads to reduction in the 
circular polarization degree of the PL. Then, we analyze 
the intricate physics of spin relaxation due to scattering 
with long- wavelength fiexural phonons and compare the 
findings with the case of graphene. Ultrafast spin relax- 
ation is predicted for free-standing and suspended SL- 
TMDs whereas for supported membranes the spin life- 
time is greatly enhanced. Finally, we discuss the relation 
between spin relaxation and charge mobility. 

Figure (TJa) shows the trigonal prismatic lattice struc- 
ture of monolayer MX2 where M (X) denotes the 
transition-metal (chalcogen) atom. Essential parts of 
the energy bands are sketched in Fig. [TJ^b). The lack of 
space-inversion center lowers the symmetry from point- 
group D^h to Dsh^ leading to spin-split energy bands and 
to complex irreducible representations (IRs) of general 
points in k-space [20]. The IRs of high-symmetry points, 
shown in Fig. [Ifb), provide useful insights to transport 




FIG. 1: (Color online) (a) Top and side views of the real 
space lattice, (b) Typical scheme of primary and satellite 
valleys in the conduction and valence bands. As implied from 
the strict one-dimensional irreducible representations of the 
double group, all bands in the K point are spin-split pO ]. 



and optical properties. We begin with the lowest-order 
processes characterized by nonvanishing phonon-induced 
scattering between valley centers, and by a temperature 
dependence that follows the phonon population (Bose- 
Einstein distribution in equilibrium) [21]. Table Ulists the 
lowest-order selection rules from which all of the nonva- 
nishing phonon symmetries are found. These modes are 
visualized in Fig. [2] for the cases of zone-center and zone- 
edge phonons. With the exception of r4 that corresponds 
to zeroth-order spin flips in the T valleys, energies of all 
other phonon modes are non- negligible ^] . Accordingly, 
charge or spin transport in the K valleys is affected by 
zeroth-order processes only at elevated temperatures. We 
summarize their properties below. 

The intravalley relaxation of either electrons or holes 
is induced by homopolar scattering due to an out-of- 
phase displacement of the two chalcogen atomic layers. 
As shown in Fig. [Sfa) and (b) for the respective cases 
of momentum and spin relaxation, the effective charge 
dipole from these long-wavelength displacements is zero 
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TABLE I: Zeroth-order selection rules in MX2 compounds. 
For spin-conserving scattering, double-group IRs are replaced 
with simpler single- group IRs. Time reversal symmetry con- 
nects K and —K points (e.g., Kq—KIq and K2—K^). 
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FIG. 3: (Color online) Direct and indirect exciton bands. The 
optical and exciton-phonon scattering processes are sketched 
under cr+ photon excitation. Real (virtual) scattering is de- 
noted by solid (dash) lines and phonon IRs are marked. Right 
blue (left red) paths lead to (j^{(T-) luminescence. 
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FIG. 2: (Color online) Atomic displacements of the 
F and K phonon modes involved in zeroth-order scattering. 



so that the relaxation is governed by a short-range scat- 
tering potential. The intravalley momentum relaxation 
is caused by thickness fluctuations of the layer due to the 
out-of-plane motion of the chalcogen atoms [Fig. [2fa)]. 
This physical picture was first identified by Fivaz and 
Mooser [14], and supported by ah initio calculations of 
Kaasbjerg et aVs who also showed a comparable contri- 
bution to the charge mobility from Frohlich interaction 
jof. Zeroth-order spin flips in the K valleys are enabled 
uniquely by homopolar in-plane optical phonons which 
do not exist in graphene structures [Fig. [2fb)]. 

The intervalley scattering between primary and satel- 
lite valleys (KxT & KxT) is relevant due to the flat 
nature of the d bands in SL-TMDs. This scattering is 
likely to facilitate the Gunn effect when applying a large 
in-plane electric field [22]. Namely, accelerated carriers 
are scattered to the satellite valleys in which the mobility 
and spin relaxation rates are different MX2 com- 

pounds with heavy (light) chalcogen atoms have a rela- 
tively small Atk (Arx) energy spacing ^], and there- 
fore can be used as n-type (p-type) Gunn diodes. 

The lowest-order processes in the one-particle picture 
are completed with the nonvanishing scattering between 
K and —K valleys. The selection rules show that such 
scattering largely affects the charge but not the spin 
transport. The lowest-order spin- flip of either electrons 



or holes is forbidden by time reversal symmetry. The 
conduction-band rule (i^axiiTI) is relevant for electron 
transport in n-type monolayers f^. The valence-band 
rule (KixKi)^ however, is likely to affect less the trans- 
port of holes due to the relatively large energy splitting 
{Ks ^ Ks or ^ Kj intervalley transitions in Fig. [1]). 

As an important application, we show that spin- 
conserving scattering between K and —K valleys is im- 
perative for understanding recent exciton PL measure- 
ments lo|-13|. Figure [3] shows the exciton bands of 
singlet composites ('bright excitons'), where upper and 
lower branches are due to the valence-band spin split- 
ting [25I. The zone-center bands comprise electron and 
hole states of the same K valley (direct excitons), where 
each of the doubly-degenerate zone-center states trans- 
form as Tq and include m/=±l excitons depending on 
light helicity. Zone-edge bands, on the other hand, belong 
to i^2,3 and comprise electrons and holes from different 
valleys (indirect excitons). In the supplemental mate- 
rial we quantify the circular polarization degree of the 
PL by modeling the absorption and relaxation processes 
that precede recombination. Here, we summarize this 
physics and show the relation with selection rules. First, 
when the exciting photon energy is between Ea and Eb 
(Fig. [3]), both upper and lower bands are excited due to 
energy broadening by impurities and substrate imperfec- 
tions [iqI, EH, 13|. Second, the phonon-assisted indirect 
absorption, a second-order process shown by dashed ar- 
rows in Fig. [3l cannot be neglected since it has many 
more available final states compared with the direct ab- 
sorption that is limited to /C = ke+k/i = Third, 
the intervalley scattering during the relaxation from the 
upper to lower F^ band flips mi without a spin-flip of the 
electron or hole [27|. The missing angular momentum is 
carried by Ki and phonons, as implied from Table |T] 
and visualized in Fig. [3] for <j+ excitation. This rapid re- 
laxation explains the measured reduction in the circular 
polarization degree [lH, and it is enabled in multivalley 
direct-gap SL-TMDs due to the (so-far ignored) unique 
coexistence of direct and indirect exciton bands. 
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In the remainder of the Letter, we study the trans- 
port due to interaction with long-wavelength acoustic 
phonons. As in other materials, the vanishing ener- 
gies of these phonons render this nonzero-order interac- 
tion important. We first mention the fundamental dis- 
tinction between charge and spin transport in monolay- 
ers that respect mirror symmetry (e.g., graphene and 
SL-TMDs). For interactions with single phonons, spin- 
conserving scattering is not affected by long-wavelength 
flexures in the out-of-plane direction [28|, whereas spin- 
flip scattering induced by flexural phonon does not vanish 
[the ZA mode in Fig. Efc)]. Mathematically, it is un- 
derstood from the scattering integral, (s/|VV|si), where 
in-plane dV /dr\\ (out-of-plane dV /dz) deformations are 
even (odd) with respect to mirror symmetry which also 
brings in ±1 for s^ = ±sj. 

To compactly express the scattering rates we employ 
the invariant method [29| . The transition amplitudes for 
spin-conserving scattering with in-plane long- wavelength 
acoustic phonons in the F, K oi T valleys are 



A,q- 



2QAvlq^ 



^A-q + S^ (^AXq)^, (1) 



where ^, A and vx are the areal mass density, area, and 
mode-dependent sound velocity, respectively, q = k/ — 
is the small wavevector difference between the initial and 
final states. Summation over tensor components is im- 
plied in the T-valley. Typical values of components of 
the deformation-potential tensor (S) are a few eV (9l.[2Qj. 
Employing the isotropic model for the phonon polariza- 
tion vectors, ^ta = {Qy^ —Qx)/q and ^la = ^./q^ one can 
resolve the contributions of in-plane transverse and lon- 
gitudinal acoustic modes. 

The spin-flip transition amplitudes from scattering 
with long- wavelength flexural phonons (^^ ~ 1) are 



' ZA,ci ' 



l kBTg{ki,kf) 
2qA 



(2) 



where 7 and 77 are mechanically-dependent parameters 
that set the flexural-phonon dispersion. Their values will 
be introduced when estimating the spin lifetimes in dif- 
ferent monolayer conditions. ^(k^,kj) can be replaced by 
Ik^+kjlg' in the zone-center, by 'Bf^q in the zone-edge, 
or by in-between, where are spin-orbit coupling 
scattering constants. The linear wavevector dependence 
in the K valley enables estimation of from spin- 
dependent energy changes in the band structure when 
applying static strain. By symmetry, shear-strain com- 
ponents of the form ej- = txz i^yz are associated with 
frozen flexural phonons, and by focusing on the spin de- 
pendent part of the static strain Hamiltonian at k = K, 



b\e+\ 



AK-&|e+|2 



(3) 



we find that a = [20| . a and b are spin-dependent 
shear deformation potentials, and is the spin split- 
ting without strain. The underlying physics is similar for 
K points in the conduction or valence bands. The value 
of needed for estimation of the spin relaxation, can 
now be readily extracted from the strain-induced spin 
splitting v'(AK + 26|e+|2)2 + |2ae+|2. Using ABINIT 
(an open-source DFT code) with Hartwigsen-Goedecker- 
Hutter pseudopotentials, we find that 0.2 eV in the 

conduction band of monolayer M0S2 |3Q|, and ~ 0.1 eV in 
its valence band [2o| . A notable feature in the conduction 
band of SL-TMDs is that can be much smaller than 
Their decoupled nature explains the difference: A^ 
is induced by small contribution from chalcogen atoms 
[13], and by interband spin-orbit coupling that in- 
volves the orbital of the transition-metal atom. 

Having values of the spin-dependent scattering con- 
stants, we quantify the i^- valley spin relaxation rate due 
to interaction with flexural phonons [Eq. ([2])]. To account 
for the spin splitting, the top valley is indexed by t and 
the lower one by £. For elastic scattering from £{t) to 
t{£)^ the spin-flip rate of the k-state is 



s-/(k) 



2m,(,)(S^^°)^fcBT 



2Fi(i,r?-l;l; 



4fcfc' 
(k+k'Y 



,(4) 



where k' = ^/mt(l^k'^ /rui^t) + (-)2m^(^) A^/n^ and m^(^) 
is the effective mass. The hypergeometric function can 
be recast to simpler forms for case-specific 77 values. The 
respective expression for T valley spin flips is similar in 
form to dl]), but with r]^r]-\-l which reflects faster spin 
relaxation (as implied from its zeroth-order selection rule; 
see Table H]). We continue the analysis and calculate the 
if- valley spin relaxation rate in two limiting cases. 

free-standing monolayers. Without a stiffening mech- 
anism to suppress violent undulations, 2D membranes 
would crumple [31]. In crystal monolayers, this mecha- 
nism is naturally provided by the coupling between bend- 
ing and stretching degrees of freedom. To lowest order 
in the coupling, the long-wavelength phonon dispersion 
is renormalized from 77 = 2to7^ = 3/2 and 7 = {/ksT j qvq 
where vq ~ Vta \/1 — v'^a/^'la [ZSI^ Hll • Assuming similar 
masses (m^ = mt = m^), the effective spin relaxation rate 
for Boltzmann distribution becomes 




Erfc(Vg^) 

1 + e-Pii 



0.1 



, (5) 



where j3K = /ksT . Due to a relatively large valence- 
band splitting, the jSx > 1 limit applies for holes at all 
practical cases where the temperature dependence is 
largely set by the complementary error function. Fur- 
thermore, in compounds with heavier transition-metal 
atoms (larger splitting), the flexural induced spin relax- 
ation of holes is slower in spite of a larger S^^. For 
example, in M0S2 where A^^ = 160 meV and in WS2 
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where =450 meV [l^, the respective spin hfetimes 
at 300 K are, respectively, ^ 0.2 ns and ^2 fis [IH. In ad- 
dition to scattering with flexural phonons, the intrinsic 
spin relaxation of holes is affected by intravalley scat- 
tering with in-plane homopolar phonons [Fig. Efb)], or 
by intervalley scattering between K and —K. Whereas 
the latter spin-flip scattering is forbidden in the zeroth- 
order by time reversal symmetry (Table H]), it is not im- 
peded by the relatively large spin splitting. Signatures of 
the homopolar and intervalley spin-flip mechanisms can 
be observed from their temperature dependence (Bose- 
Einstein distribution of the involved phonons). 

The spin relaxation of electrons is much faster due to 
the small spin splitting in the conduction band. In M0S2 
where A^^^ = 4 meV [Tol . the room-temperature spin life- 
time is Ts ~0.05 ps [32|, and it increases noticeably only 
below 50 K. Interestingly, the spin relaxation of electrons 
is enhanced in compounds with lighter chalcogen atoms 
(their weaker spin-orbit coupling leads to smaller A^^). 
Furthermore, the spin relaxation diverges in the patho- 
logical limit of A^^'' [33| . Finally, repeating the anal- 
ysis for suspended monolayers by employing r] = l and 
7 = vo^/i, we find that < 1 ps at room-temperature 
even for suspension-induced strain levels of e ~ 1%. 

supported monolayers. Another means to stiffen the 
membrane is naturally provided by van der Waals (vdW) 
interactions when the monolayer is placed on a substrate 
[3^ . The support brings in a minimum cutoff energy 
for out-of-plane displacements. The cutoff energy, Qc 
= h y^/^s/m-^, is calculated from the average vdW inter- 
atomic force constant between the monolayer and the 
substrate (a^s), and the average atomic mass of the mono- 
layer {Mr,). In the long- wavelength limit, we can therefore 
approximate the dispersion of flexural phonons by 77 = 
and 7 = l^c/^- The effective spin relaxation rate becomes 



with flexural phonons follows (|2]) where 



1 _ (2m^S-)' 4 + 2pK 



(l + efe)/32 



(6) 



Pc = ^c/kBT <l is assumed. The temperature depen- 
dence is quadratic for f3c < Pk < 1 and exponential for 
/3c < 1 < The substrate coupling brings in slower re- 
laxation which at room temperature reaches Ts^T ps 
(~ 1 ns) for electrons (holes) in supported M0S2 with 
= 1 meV [35] . The spin-lifetime enhancement from 
?7 = 3/2to7^ = 0is sharper for electrons due to the small- 
ness of A^^ compared with A^\ As a result, phonons of 
longer wavelength are involved in energy-conserving spin 
flips in the conduction band. 

We compare the spin relaxation induced by flexural 
phonons in graphene and MX2. The space inversion 
symmetry in graphene mandates spin-degenerate energy 
bands, resulting in anisotropic relaxation that depends 
on the spin orientation. Deflning the latter by out-of- 
plane polar angle (^s) and in-plane azimuthal angle ((/)s), 
the spin-flip transition amplitude due to elastic scattering 



-cos2^e^^+ 

- 1 kx 



•(7) 



The other angles, (/)± = ^(tan -"^-^^itan ^ -^)^ are 
due to the Dirac-cone energy dispersion. Formal ana- 
lytical derivation of this result will be presented in a fu- 
ture long publication. Here we mention that the pref- 
actor k s\v? (j) = q'^ / , previously found by numerical 
techniques [36j, originates from space-inversion symme- 
try and spin-dependent energy dispersion in the cone re- 
gion. Away from the Dirac point (/c=0), the spin re- 
laxation timescales are longer in graphene than in SL- 
TMDs on accounts of the higher power-law in phonon 
wavevector {q^ vs q) and the relative smallness of in 
carbon-based systems (~ 10 meV [37|). 

All these findings lead to a counterintuitive relation 
between charge mobility and spin relaxation in 2D mem- 
branes. Whereas increased stiffness has been shown to 
be associated with slower spin relaxation, its coupling 
with charge mobility seems to have the opposite trend. 
For example, high mobility in supported membranes is 
a token of diminished effect from adsorbents and sub- 
strate imperfections. A smaller coupling of the mem- 
brane to such parasitics would enable freer and softer 
out-of-plane undulations. Therefore, an inverse trend 
between spin and momentum relaxation, a hallmark of 
Dyakonov-Perel spin dephasing processes, can be real- 
ized in a Elliott- Yafet spin flip system (scattering with 
flexural phonons). This physics is universal in relatively 
clean 2D membranes that respect mirror symmetry. In 
such membranes, charge transport is decoupled from har- 
monic out-of-plane undulations while spin relaxation is 
not severely affected by the presence of impurities. 

In conclusion, we have presented a transparent theory 
of intrinsic transport properties in single-layer transition- 
metal dichalcogenides. Lowest-order scattering processes 
were identified in both charge and spin transport regimes, 
and were found relevant to the transport of free-carriers 
at elevated temperatures and for the energy relaxation 
of hot excitons. As a group- theoretical study, findings 
of previous important theories were mentioned as special 
cases. In addition, the spin relaxation induced by scat- 
tering with long- wavelength fiexural phonons was quanti- 
fied. For electrons, the ultrafast rate is attributed to the 
typical softness of 2D membranes and to the small spin 
splitting in the conduction band. The relatively large 
splitting in the valence band, on the other hand, renders 
p-type monolayers better candidates for preserving spin 
information. 
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